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Abstract 

Extended Einstein-Maxwell model and its application to the problem of 
evolution of magnetized Bianchi-I Universe are considered. The evolution of 
medium magnetization is governed by a relaxation type extended constitutive 
equation. The series of exact solutions to the extended master equations is 
obtained and discussed. The anisotropic expansion of the Bianchi-I Universe is 
shown to become non-monotonic (accelerated/decelerated) in both principal 
directions (along the magnetic field and orthogonal to it). A specific type 
of expansion, the so-called evolution with hidden magnetic field, is shown to 
appear when the magnetization effectively screens the magnetic field and the 
latter disappears from the equations for gravitational field. 
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1 Introduction 



Observational cosmology and astrophysics indicate that many subsystems of the 
Universe, which have different length-scales: planets, stars, galaxies, clusters and 
superclusters of galaxies, possess inherent magnetic field (see, e.g., the reviews [Il[2|, 
[3] and references therein). The magnetic field interacts with material environment 
and this results in the magnetization of the latter |3]. According to basic Einstein's 
ideas, the stress and energy of all three constituents: medium, magnetic field and 
magnetization, - act as sources of gravity field. It seems to be reasonable from 
physical point of view, that for majority of astrophysical and cosmological objects 
the contribution of a material medium, as a source of gravity field, dominates over 
the contributions of magnetic field and magnetization, nevertheless, the stress and 
energy of magnetization itself can be comparable with those of pure magnetic field. 
In other words, when we deal with exactly integrable models, for which the magnetic 
field is considered as a non-negligible source of gravity field in comparison with 
the contribution of the material medium, one should also take into account the 
magnetization as a cross-effect. 

Cosmological Bianchi-I model seems to be the most convenient for testing this 
idea due to its three specific features. First, this model is non-stationary, all three 
constituents, medium, magnetic field and magnetization, can evolve with different 
rates, thus displaying their dynamic non-equivalence. Second, the Bianchi-I model 
belongs to the class of the exactly integrable ones, thus providing the analysis of sin- 
gular behaviour of the model. Third, the Bianchi-I model is spatially anisotropic and 
thus admits a self-consistent description of the uniaxial configuration of the dynamic 
system containing interacting matter, magnetic field and magnetization. There is a 
number of anisotropic cosmological models, in which magnetic field is accompanied 
by perfect fluid (see, e.g., [5], [6] and references therein), or by non-equilibrium (vis- 
cous) cosmic fluid (see, e.g., [3 El [9]). If in addition to matter and magnetic field 
the third "player" , namely the magnetization, appears in a cosmological dynamics, 
one can expect that the rate of evolution of the Universe modifies. Our expectations 
are motivated by the analogy with dissipative phenomena in cosmology, described 
in the framework of causal (extended) thermodynamics [10] - [21]. In this theory 
the extended constitutive law for the viscous fluid contains a time parameter, which 
is known to introduce a specific time scale into expansion rate. Since the evolution 
of magnetization has usually a character of relaxation, a new time parameter, say 
relaxation time, r, also must appear and establish a new time scale. The inter- 
play between r and expansion rate parameter(s) can introduce a qualitatively new 
aspects into cosmological dynamics. 

In order to formulate a self-consistent Einstein-Maxwell model, taking into ac- 
count the polarization and magnetization of a non-stationary material medium, two 
key elements are necessary. The first one is an adequate energy-momentum ten- 
sor of the electromagnetically active medium, which forms a source term in the 
right-hand-side of the Einstein equations. Here we derive explicitly such a tensor 
using Lagrangian formalism for the stationary non-conducting medium with uniax- 
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ial symmetry in case when magneto-electric cross effects are absent. Expressed in 
terms of induction tensor and Maxwell tensor, this quantity happens to be a sym- 
metrized Minkowski stress-energy tensor (see, e.g., |25]-[29] for a review, historical 
details and terminology). When a medium is non- stationary, the required effective 
stress-energy tensor is assumed to have the same formal structure. Nevertheless, the 
induction tensor acquires now a new sense: it is a sum of the Maxwell tensor and of 
a polarization-magnetization tensor, the latter quantity being considered as a new 
dynamic variable. The second key element of non-stationary model is the extended 
constitutive equations, which establish relations between polarization-magnetization 
of a medium and electromagnetic field strength. We use here the simplest extended 
constitutive equations of a relaxation type, which are formulated phenomenologi- 
cally, based on the well-known analogs from causal (extended) thermodynamics [10] 
- [23]. The main goal of this paper is an application of the formulated extended 
Einstein-Maxwell (EEM) model to the description of the magnetization dynamics 
in the Universe, considering the corresponding master equations in the context of 
Bianchi-I anisotropic cosmological model. 

The paper is organized as follows. In Section 2 we briefly discuss the principal 
details of the EEM- model. Particularly, we introduce the effective stress-energy 
tensor describing the electromagnetic field and the polarization - magnetization in a 
medium, as well as we introduce the constitutive equations of a relaxation type. In 
Section 3 we adopt the EEM model for the symmetry related to Bianchi-I cosmolog- 
ical model and discuss the reduced Maxwell, Einstein and constitutive equations. In 
Section 4 we consider exact solutions of the dynamic equation for the magnetization 
with variable relaxation time parameter, discuss general properties of these solu- 
tions and some interesting particular cases. In Section 5 we obtain exact solutions 
of the Einstein equations for the case of variable relaxation time parameter. We 
distinguish two principal submodels in this context. The first submodel describes 
a paramagnetic/diamagnetic dust (Subsection 5.1) and contains three exactly inte- 
grable particular cases. The second submodel describes the so-called longitudinal 
quasi- vacuum (Subsection 5.2) and also contains three exactly integrable particular 
cases. In Section 6 we obtain exact solutions of the EEM- model for the case of 
constant relaxation parameter. In Subsection 6.1 we establish a law of the magneti- 
zation evolution. In Subsections 6.2 we consider the model with hidden induction. In 
Subsection 6.3 three particular cases of exact solutions of the EEM model, describing 
the submodel with vanishing total longitudinal pressure, are discussed. In Subsec- 
tion 6.4 the example of cosmological dynamics with non-homogeneous non-linear 
equation of state of the magnetized matter is studied. Discussions form Section 
7. Appendix contains the description of the procedure of variation of the tetrad 
four- vectors, which is used in Section 2. 
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2 Extended Einstein-Maxwell model 



2.1 Stationary media with uniaxial spatial symmetry 

Let us consider a preliminary model with the following action functional 

•S" = y d^Xy/^ I h -^v(matter) + '^C'^"''^ FikFmn^ , (1) 

as a /imt for construction of the extended Einstein-Maxwell model. Here R is the 
Ricci tensor, g is the determinant of Qik, A is cosmological constant, k, = StcG, 
G is the gravitational constant. The quantity Fmn is the Maxwell tensor, Fmn = 
dmAn — dnAm, Am is a potential four-vector of the electromagnetic field. L(matter) 
is the pure medium contribution to the Lagrangian, we assume that this scalar 
does not depend on Fmn- The quantity (7*'^™-" is the linear response tensor, which 
describes the influence of matter on the electromagnetic field. This tensor is assumed 
to possess the following symmetries 

^ikmn ^kimn ^iknm ^mnik 

and is a function of the metric gmn-, time-like velocity four-vector of the medium as 
a whole, , and some space-like vector , pointing the privilege direction in the 
medium. These four- vectors are orthogonal and normalized by unity, i.e., 

g„WX^ = Q, g„UPU'' = l, g„X^X'^ = ~l. (3) 

When C^^™"^ contains the Riemann tensor, the Ricci tensor and the Ricci scalar. 



we deal with non-minimal Einstein - Maxwell theory [SUl EI]- When C^^"^^ in- 
cludes the covariant derivative of the velocity four- vector ViUk, the corresponding 
model describes dynamo-optical effects [32]. It is worth stressing that here we re- 
strict ourselves by the case when C'**^'"" = C*'^™'"[(y'pg, X']. The procedure of 
phenomenological reconstruction of the material tensor C^^^^ is well-known [H [33] . 
In order to obtain C**'^™" as a function of gmm and X^ for the medium with 
uniaxial symmetry one uses, first, the standard decomposition 

/^ikmn ^TTQTT^ i s:ik jrmn Ip ik mn (,,—l\lp i , ,lp ^ik jrmn i mn jrik \ ( a\ 

— ^ \%^«s ^ ~^ pq^ ls{f^ ) +^ ls"pq +^ ls"pq\j ' 



where dielectric permittivity tensor, e^, magnetic impermeability tensor, (/i 
and magneto-electric tensor u''^ are used, defined as 



The quantity 5*^ is the Kronecker tensor, e*'^-'^ is the Levi-Civita tensor 
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j^ikjs jg ^YiQ completely skew - symmetric Levi-Civita symbol with i^o^^s _ r^^ie 
Levi-Civita tensor provides the dualization procedure: F*i^ = ^eikmnF"^"' . The sec- 
ond step of reconstruction of (7*'^™" tensor for the case of uniaxial symmetry is a 
phenomenological representation of the tensors e''^, and u''^. In the simplest 

case, when the magnetoelectric cross-terms are absent (i/'^ = 0), such relations are 

e^P = e±A^P+(e^-e\i)X^XP, (/x-^)'^' = — A'^ + | — - — ) X'X^ . (7) 

Here the scalar quantities e± and are transversal and longitudinal coefficients 
of dielectric permittivity, respectively, scalars fi± and fi\\ represent coefficients of 
transversal and longitudinal magnetic permeability, respectively, A'^ = g''^ — U^U^ 
is a projector. Finally, (jl]) can be rewritten as 



C 



ikmn 



1 



1, 



g^mgkn _ ^^n ^km^ ^ ^^^^^ _ ^) ^k^ ^mn ^,1 ^J, 



+ 



2 



( ^ ^ \ Siik xmn ( ^ ^ | Ak mn 

[e^ - e\\] di, -[ — - — j^ Is 



(8) 



When e± = e\\ = e and fi± = fi\\ = fi, the obtained tensor of material coefficients 
(jikmn QQYgrg well-kuowu isotropic one (see, e.g., 133]). 

2.1.1 Maxwell equations 

Variation of the action functional ([1]) with respect to the four-vector of electro- 
magnetic potential Ai gives the Maxwell equations with vanishing current of free 
charges 

Vfc(C*'=™"F^„) =0. (9) 
In this case the induction tensor, if*^, is equal to the expression in the parentheses. 

2.1.2 Equations for gravity field 

Variation of the action functional ([1]) with respect to metric gik yields 

JDik 1 Ak p A ik , ,.T^«fc rpik rpik , rpik /-i ri\ 

-n- — 2^ ^ ^ ^ (total) 1 ^ (total) — ^ (matter) + ^ (eff) • \^^) 

The symmetric stress - energy tensor of the material medium T^'^atter)' defined as 

^(matter) = 7^77" ("s/^-^' (matter)) , (H) 
V 9 ^9ik 

can be written in the standard form 

^(matter) = WWU'' + q'U^ + q'^U' - PA''' + rf^ , (12) 
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where W is an energy density scalar of the matter, is a heat-flux four-vector, 
P is the Pascal pressure and 11*^ is an anisotropic pressure tensor. The effective 
stress-energy tensor has the form 

Tl^) = -7A=Fp,F^n^ (v/=^C") , (13) 

it can be calculated directly for the medium with uniaxial symmetry using the 
decomposition ([8]) and the following formulas for the variations 5U^ and dX'^: 

5U' = ^6g.k {U'g^' + U'g") , 5X^ = ^6g,k {x'g'' + X'g'') . (14) 

The grounds of the formulas f|T^ are presented in the Appendix. The variation 
procedure yields 

4ff) = \9"'C^'"'''Fp,Fmn - 1{C'"''"F'^ + C'"'P'^F'JF„ , (15) 
i.e., the effective stress-energy tensor is explicitly symmetric and traceless. 



2.1.3 Resume 

When the non-conducting medium with uni-axial symmetry is the stationary one, 
the Einstein-Maxwell model consists of three ingredients: 
(i) Maxwell equations 



ik 



0, VkF 



*ik 







(ii) constitutive equations 



(16) 



H 



ik 



c 



ikmn ; 



(Hi) gravity field equations 



1 . 



g'^R = Ag'" + k WWW' + q'W + q'^W - PA''' + U 



.Ak 



ijrk I „irTk 



-k T ri 



, ik 



lik 



+ 



+ K 



1 

_^ikTT rpr. 
^g ^mn^ 



-{W^F'^ + H'^F\) 



jkm ; 



(17) 



(18) 



It is the exact result, which follows from the variation procedure. The electromag- 
netic part of the total stress-energy tensor in the presented form coincides with the 
symmetrized Minkowski electromagnetic energy- momentum tensor (see, e.g., the re- 
view [29j). It is manifestly symmetric, traceless and does not depend explicitly on 
the choice of W. 
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2.1.4 Extended Einstein-Maxwell model 



When the medium is non-stationary, the Einstein-Maxwell model should be mod- 
ified accordingly, the stationary equations (fT6l)-(fT8l) can be used as a hint. Our 
ansatz for such extension is the following: we consider the electrodynamic equations 
in the same form (fT6!) . but express the induction tensor in terms of polarization- 
magnetization tensor M*'^; instead of stationary constitutive equations (ITTl) we 
introduce extended constitutive equation of a relaxation type for the tensor of 
polarization- magnetization; the gravity field equations have the same form (fTSll . 
but W'' is replaced by F*'^ + M^^, i.e, the basic set of master equations for the 
Extended Einstein-Maxwell model is 

Vfc [f''' + M'^) = , VkF*''' = , (19) 

tDM''' + M'^ = x'^'^'^Fmn , (20) 

+ 



jk I „iTTk I „kTTi u \ik , Ttik 



^ ifc P rpmn rpim rpk 



+ K 



(21) 



Here r is a relaxation time and D = U'^Vk is a convective derivative, is a 

linear susceptibility tensor, which can be expressed in terms of C***^™" as 



ikmn 



^ (jikmn _ 1 ^ginigkn _ ^in gkm\^ _ ^^2) 



The quantity M'^ is considered in this context as a new variable. Such extension of 
the constitutive equations can be regarded as phenomenologically motivated, when 
the well-known analogs from rheology and extended thermodynamics pLQj - [24] are 
taken into account. Thus, the novelty of the presented model is connected, first, with 
the extended constitutive equation fl20l) . second, with the modified electromagnetic 
source in the gravity field equations (see the last line in (l2T])). 

3 Bianchi-I Universe 

Consider the rotationally isotropic Bianchi-1 cosmological model [5], [6] with the line 
element 

ds^ = de - a\t) [{dx^ + {dx^] - c^{t) {dx^f . (23) 

We suppose that the velocity four- vector of a matter, U^, has the form = 6^ 
and the space-like vector X* is X* = 5\/c{t). They satisfy the relations DU' = 
and DX' = 0. The symmetry of the model prescribes the non-diagonal components 
of the total stress-energy tensor to vanish. In this sense to make the model self- 
consistent one can consider, for instance, the case when the magnetic field four- 
vector and the magnetization four-vector are parallel to the x^ axis. For such a 
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configuration of tlie electromagnetic field the total stress-energy tensor, T^^otai)' 
four non-vanisliing components: 



^oVtal) = = W + X , T^\total) = ^2Vtal) = "^(tr) = "(^(tr) + X) , 

7f(totai) = -Vii = -{Pii -X), X^ \h'^F,2 = \[F'^F,2 + Mi^Fi2] . (24) 

Here P(tr) = P{i) = P(2) and P\\ = P(3) are tlie transversal and longitudinal pres- 
sure components, respectively, coinciding with the corresponding eigenvalues of the 
material part of the total stress-energy tensor. 

3.1 Reduced Maxwell equations 

We assume that the solutions inherit the spacetime symmetry and the quantities Fj^ 
and Mjfc are the function of cosmological time only. Thus, it follows from the second 
subsystem of Maxwell equations f|T6|) that the spatial components of the Maxwell 
tensor are constant, i.e., F^p = const. The only non-vanishing term is F12 = const, 
since the magnetic field points along the axis. Note that the tetrad component 
of the magnetic field, B(t), is connected with the constant Fu by the following 
relationship 

We assume, that the electric field, electric polarization and magneto-electric cross- 
effect are absent and the medium is locally neutral. These requirements guarantee 
that the Maxwell equations are satisfied identically. Thus, in the model under 
consideration the evolution of the magnetization is governed only by the constitutive 
equation (l20l) . 

3.2 Reduced constitutive equations 

When W = 6q in the metric ( l23l) the acceleration vector DU^ vanishes, yielding 
jjjj^ki _ £)(^g«'='*f/^) = 0. The magnetization four-vector has only one component: 

M' = M*'^Uk = -X'M , M''' = -7]'^^ Ml = r]'^^XiM . (26) 

Extended constitutive equations (120!) reduce to one equation of the relaxation type 

rM(t) + M=fl-l)i?o^, (27) 
where the dot denotes derivative with respect to time. 



8 



3.3 Reduced Einstein's equations 

The gravity field equations reduce to the following system 

2 



a\ a c 
- +2— =A + k{W + X) 
a J ac 



a c 



- + - + —= A- K{P^tr)+X) 



a c 



a \a/ 

Differentiation of Einstein's equations leads to the conservation law 



W + 2{-] W + P( 



(tr) 



^ C 



W + P\\ + = 0. 

a, 



(28) 

(29) 
(30) 

(31) 



Note that X and its derivative enter the conservation law with the multiplier, which 
does not depend on c(t) (on a{t) only). The function X(t) is now 



Xit) = -B{t) [M(t)+i?(t)] = ;i?, 



2 (^'^{to) 
° a4(t) 



1 + 



a\t) M{t) 
a^{to) Bo 



(32) 



When the magnetization M vanishes, the quantity X(t) is non-negative. Never- 
theless, when M is non- vanishing, X{t) can be negative during some time interval, 
or be equal to zero. The total energy density W = W + X is assumed to remain 
positive. 



4 Magnetic relaxation with variable time param- 
eter T{t) 

The relaxation equation (!27|) and the quantity X{t), the source of the gravitational 
field (1321) . contain the function a{t) and do not contain c{t). This means that it 
is reasonable to split the master equations into two subsystems, dealing with the 
longitudinal (P||, B{t), M{t)) and transversal (-P(tr)) quantities, as well as, with 
functions c(t) and a{t), describing the evolution of the Universe in the longitudinal 
and transverse directions, respectively. 

Notice that the sign of the right-hand-side of the equation fl27|) depends on the 
sign of the expression (l//i|| — 1). According to the standard definitions, when > 1 
we deal with the so-called paramagnetic medium; in this case the expression in the 
parentheses in ( 1271) is negative. When < < 1, the medium can be indicated as 
diamagnetic, this case relates to the positive expression in the parentheses in fl27j) . 
Finally, when >> 1, one can say that the medium is in a ferromagnetic phase. 
We do not consider the cases with negative and vanishing In principle, can be 
treated as a function of time. This means that the sign of the expression (l//i|| — 1) 
may change with time. This problem itself is very interesting, but the numerical 
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calculations are needed to solve the corresponding master equations. We focus here 
on the analytical solutions, and do not consider the case of variable 

Generally, the master equations of the EEM-model are self-consistent. In order 
to find M{t) from f l271) we have to know a{t). In order to find a{t) from (l30ll 
we should know A'(t) which includes M{t). To solve such a problem one should 
introduce some extra ansatz. We prefer to start from the equation fl27|) . and our 
ansatz concerns the function T{t). There are at least three ways to introduce the 
relaxation time parameter T(t). The first approach (which is the simplest one) 
assumes that r is constant. We consider the model with r(t) = tq in Section 6. 
In the second approach the relaxation time T{t), the bulk viscosity coefficient C(^)) 
etc., are modeled as a power-law functions of the energy density scalar W (see, e.g., 
p!7] - [21]). We do not consider such options in this paper. In the third approach 
r(t) is assumed proportional to the inverse rate of expansion. For instance, in the 
isotropic Friedmann model {a{t) = c{t)) the new dimensionless variable ^ = TH{t) 
is frequently used, where H(t) = a/a is the Hubble expansion parameter (see, e.g., 
|34] ) . In this section we consider the third version of the representation of T{t). 

Since only the function a{t) enters the relaxation equation and the function X{t), 
we consider r(t) = . It is convenient to use a new variable x = and a function 



a{to) 



I (a) 



X 



-, for which 



- = H, 



(a) 



X] 



+ ^(a) 



The three-parameter family of solutions to (I27|) reads 

^0 (I 
(1 - 20 Ull 



M(x,e,M(to),/i||) = Mox" 



« + 



1 



X — X 



Here 2^ 7^ 1, the special case ^ 
corresponding expression for X {x,^, Mq, is 

1 



I will be considered in Subsection 4.2. 



X 



where 



2^/^11 



20 



and ir2(e,Mo,/i||) + 



(/^ll-l) 
/i||(l-20 



(33) 



(34) 



The 



(35) 



(36) 



are constant. We use for simplicity the definitions Mq = M(to) and Bq = B[to), and 
below we omit the parameters ^, Mq and fi\\ in the arguments of M, X and Ki,K2. 



4.1 Non-resonant magnetic relaxation: ^ > 0, ^ 7^ | 

For such ^ the magnetization decreases when the Universe expands, i.e., M{t 
00) — >■ when a{t) increases. The relaxation of the initial magnetization Mq is 
characterised by the function x~^^^. The relaxation of the induced magnetization 
Fi2(l//U|| — 1) is characterised in the limit x ^ 00 by the function x^^^^, when ^ > |, 
and by the function x"^, when ^ < ^. 
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4.1.1 Extremums and zeros of the function X{t) 



The formula ( !35l) shows, that the behaviour of X{x) is monotonic when K1K2 > 0. 
When K1K2 < 0, there exists one extremum at t = such that 




A',(2^ + l)' 2f,||(2? 

In addition, the inequality t^, > is valid for 4^ li^Til > {2^ + 1) \K2\. It is a require- 
ment to the ratio Mq/Bq, and we suppose it is valid. The extremum is a minimum 
with negative value X{t^), when | < 2/i|| and a maximum with positive value X(t^), 
when I > 2/i||, for both paramagnetic and diamagnetic medium. Taking into ac- 
count the definitions of Ki and K2, fl36|) . one can find the following possibilities for 
X[t) to have an extremum. When the medium is paramagnetic, i.e., > 1, one 
obtains three different situations. 

(i) I < 2, < 0. 

In this case there is a minimum. The curve X{t) passes its zero- value point, if 
X{to) > 0, and tends asymptotically to zero as X{t) cx a^^~« , when a{t) 00. 

(ii) I > 2/i||, K2 < 0. 

This case corresponds to a maximum. The curve X{t) passes its zero- value point, if 
X{to) < 0, and tends asymptotically to zero as X(t) oc a~*^, when a(t) — > 00. 
(m) 2 < i < 2/i||, i\:2 > 0. 
This case is analogous to the first one. 

Likewise, for the diamagnetic medium < 1 one obtains the following situations, 
(z) I > 2, K2 < (maximum). 
(ii) I < 2/i||, K2 < (minimum). 



Hi) 2 > I > 2/i||, > 



maximum 



Besides, there exists a time moment, t^^,, when the source X vanishes. According to 
( 153]) the condition A'(t,,*) = can be satisfied, when 




but this is possible if the constants Ki and K2 have opposite signs. Thus, the 
existence of zeros of A'(t) assumes the same requirements K1K2 < 0, as conditions 
for the existence of extremums. 



4.1.2 Asymptotic case: fast relaxation, ^ 

When ^ ^ the formulas (IMI) and (153]) give 



-2/1 i\ -v/^N ^^(^o)__4 



M(x)|5_o = 5o x-M 1 , A" (x) 15^0 = ^7^^"' • (39) 

In a paramagnetic medium the ratio M{x)\^^q/ Bq is negative, in a diamagnetic 
medium it is positive. In both cases X{x)\^->o is positive. 
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4.1.3 Asymptotic case: slow relaxation, ,^^00 
When ^ — i> cx), one obtains 

M(x)|^^oo = Mo , A'(x)|5^oo = [Bo + Mox^) . (40) 

In these formulas fi\\ is absent. 

4.1.4 Particular case: 2^/i|[ = 1 

Studying the extremums of the function A'(t) we have found that the condition 
2^/i|l = 1 (or, equivalently, Ki = 0) gives some threshold value for the parameter ^. 
Let us consider this particular case separately. When 2C,fJ^\\ = 1 

M{x)u^^ = (Mo + So)x-'^ii - Box-^ , (41) 

2n| 

and the function X{t) monotonically tends to zero as x 00 

X{x\^=^ = \b,{B, + Mo) a;-2(^+''ii) ^ X, x<'^^^\) , (42) 

remaining positive or negative depending on the sign of its initial value Xq. Obvi- 
ously, the value X{x)^^^i^^^ is positive or negative depending on the sign of the sum 
Mo + -Bq) and is equal to zero identically when M(to) = —Bito). In the last case 
the information about magnetic field and magnetization does not enter the Einstein 
field equations and we deal with the so-called "hidden" magnetic field [31]. Notice 
that in the ferromagnetic phase ^\\ » 1, which corresponds here to the model of 
fast relaxation ,^ << 1, the source X{t) disappears very quickly. 



4.1.5 Particular case: vanishing initial magnetization, Mo = 



This particular case can be studied in more detail. The minimum of the function 
X{t) exists at the point given by 



4e(2e/i| 



11 



.«(^o). 



1 



(43) 



(2e+l)(/"|| 

Since we put Mq = 0, we should check especially the condition a{t^.) > a(to). When 
the medium is paramagnetic, this condition is satisfied for ^ > |. This corresponds 
to the first case of the classification of the extremums: the function X{t) starts from 
the positive initial value Xq, takes its zero value at t**, reaches its minimum at 
and tends to zero asymptotically, remaining negative. Here 




mm - 1) 
im - 1) 



> 1 



2^ll(2e+l) [(^||-l)(2e + l) 



4€ 

i-2e 



< 0. 



(44) 



(45) 



The inequalities for the diamagnetic medium can be obtained analogously. 
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4.2 Resonance magnetic relaxation: ^ 

When ^ = |, the special solution to (|27|1 exists 

' 1 



M{x) 



X 



Mn + 2Bo 



1 los,x 



(46) 



which tends to zero at t ^ oo (x ^ oo). The term "resonance" relates to special 
case when the relaxation parameter r(t) coincides with half of the characteristic 
expansion time l/if(a), where if(a) = ^. The function X{x)^^^i reads 



X(x) 



-Bn X 



B0 + M0 + 2Bo 



1 loffx 



(47) 



The function X(x] 



reaches its extremum value 



B\U) ( 1 



when 



log 



'a(t* 



^a(to) 



i + 2(;.„-l) 



Mo 

-Do 



(48) 



(49) 



In a paramagnetic medium this extremum is the minimum, the condition a(t*) > 
a(to) assumes that Mq/Bq > (1 — 3fi\\)/2fi\\. When the medium is diamagnetic, 
the extremum is the maximum, the condition a(t*) > a(to) is satisfied if Mq/Bq < 
(l-3/i||)/2^l|. 



4.3 Magnetic instability: ^ < and Mq 7^ 



In this case the magnetization 



Ba 



X — x^ 



(50) 



increases as X"^/'^'. In the special case Mq = ^"(''14,^^1') the term disappears, and 



the magnetization decreases as 2; ^. The function 

1 



i^'i + K2 x~"^^ 



(51) 



decreases when ||| < 2 and increases when yij > 2. When yij = 2 it tends asymptot- 
ically to the constant value 



X (cx)) = -Bo 



Mo + Bo 



2/ii| 



(52) 
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5 Cosmological evolution in case of variable re- 
laxation time 



Below we discuss exactly solvable models of four selected types, describing cosmo- 
logical evolution of the magnetizable medium. In order to explain our choice let us 
note that when X depends on time via a{t) only, the equation (!30|) does not contain 
the information about c(t), if P\\ also depends on a only, i.e., P\\ = P\\{a{t)). The 
vacuum- type equation of state W + P\\ = belongs to this class of models, since c{t) 
disappears from conservation law (j3Ti) . and its solution of the form W{a) = —P\\{a) 
pics out fl30|) from other Einstein's field equations. 



5.1 First example of cosmological evolution: 
paramagnetic / diamagnetic dust 

When the medium behaves as a dust, i.e., P\\ = P(tr) = 0, the optimal way to solve 
the master equations is the following. First, we solve (|30l) . which in terms of variable 
X takes the form of the equation for H?s (x) 



d 



X 



Tx (^^^i 



A + -kBI 



Ki x-^ + K2 a;"(^+t) 



(53) 



Second, we search for the function a(t) using the quadratures 

dx 



± 



"(to) 



(54) 



When a{t) is found, it is convenient to search for c{t) using the substitution Y{t) = 
c(^\ ali^- equation governing the evolution of Y{t) follows from (!29|) and (!30|l 



Y{t) + Yit) 



-h + ^KXiait)) 







(55) 



and is a linear differential equation of the second order with coefficient depending 
on time. The initial data for Y{t) follow from the definition and from the equation 

m 

A + k{Wo + Xo) 



Y{k 



Yit, 



(56) 



Finally, we search for W{t) using fl28|) with obtained a{t), c{t), a and c. In the 
process of integration of the equation fl53l) a new resonance value of the parameter 
^, namely ^ = 1, appears. Indeed, when ^ 7^ 1, the solution is 



(a) 



X] 



A 
^3 



- + K^x' 



Ki x-^ + K2- 







X 



(57) 



14 



where 



A kB^ 



Ki + K2- 



1-0 



(58) 



When ^ = 1, i.e., the relaxation time parameter r(t) coincides with ^, we should 
replace fl57j) by 



where 



/A ~ 



(59) 
(60) 



5.1.1 General properties of solution (t) 



Asymptotics 

The quantity plays a role of Hubble function, describing the Universe evolution 
in the plane x^Ox^. When A 7^ and | > —2, the (positive) asymptotic value of 



this function, Ha( oc) 



3 ' 



is well-known for the Friedmann isotropic model. The 



corresponding asymptotic behaviour of a{t) is a{t) oc exp{y|- t}. When A = 0, the 
asymptotic formula for H(^a) is predetermined by the value of the parameter ^: 

a) if I < 1, if(a) oc x'^^^^^ and a{t) oc 

b) if I > 1, H(^a) oc x~^, and a{t) oc ti. 

When I = — 2 the last term in the expression for H(^a) does not depend on time 
and redefines effectively the cosmological constant 



A A* 



Ah 



Mo ^ - 
Bo 2 fill 



1 



(61) 



When I < —2 the last term in the expression for H(^a) fl57|) is the leading order term 
at if: — 00, and the Universe collapses in the cross-section x^Ox'^. It is an exotic 
case, and we do not focus on it. 



Extremums 

Generally, the behaviour of the H(^a){x) function is non-monotonic, i.e., there are 
intervals with if(a)(t) < as well as with H(^a){t) > 0. The necessary condition for 
the existence of extremums marks the points X(i),X(2), in which 



K?. X(^s) - 



6 



2- 



2^+1 

(i-0"(^) 



0, 



(62) 



Formally speaking, this equation can give (m) real roots, the function H(^a) can 
possess (m) extremums, and, consequently, (m) points can appear, in which the 
transition from the accelerated expansion to the decelerated one (and vice-versa) 
takes place. To illustrate these possibilities, consider now several particular cases. 
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5.1.2 Ki = K3 = 

Let us choose the initial values if(a)(to), Mq, F^, a(to), as well as, ^ and fi\\ param- 
eters in an appropriate manner to provide the relations Ki = and = 0. It is 
possible when 



For this case X(x) is the monotonic function 



Bn 



(63) 



(64) 



and H(^a){x) simplifies essentially 



X 



-2(1+M||) 



(65) 



Here 7^ | to avoid the relation ^ = 1. For such H(^a){x) the equation (|54l) can be 
easily integrated 



a{t) 
a(to) 



cosh 



;i+/iii)(t-to) 



^(a)(io)\/^sinh 



:(l+/i||)(t-to) 



1 + M|| 



Hu){t) takes an explicit form 



Hia){t) = 

The asymptotic behaviour of this solution is given by 
a{t 00) =^ a(to) 



H{a)ito) + y f 


tanh 




h/i||)(t-to) 


+ if(a)(to) 


tanh 




h/i||)(t-to) 



(66) 



(67) 



Ml + i^(«)(to)J^ 



l + MI 



exp J^(t-to) 



H(^a){t 00) 



(68) 
(69) 



Returning to the function c{t), note that the equation fl55|) with a(t), given by fl66l) . 
is the Hill equation [35j with imaginary argument. Since X{t —>■ 00) 0, the 
appropriate asymptotics of Y{t) and c{t) are 

y (t ^ 00) oc exp <{ ^ 1 1 , c(t ^ 00) oc exp I W ^ t ^ , (70) 



thus, the Universe isotropizes. To specify c{t) consider three special examples. 
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First special example: if(a)(to) = 

When if(a)(to) = \f^, X = 0, Mq + Bq = and the behaviour of a{t) is governed by 
the de Sitter expansion law 



a{t) = a{to) exp { \l ^{t -to)} . 



(71) 



The solution for c{t) is 
kWo 



c{t) 



c{to) 



2Aa2(to 



exp 



'-{t-to) 



2Aa2(to 



exp 



it -to) 



(72) 

The asymptotic behaviour of c{t) at t — oo is the same as for a{t). This solution 
corresponds to the model of hidden induction, mentioned in the subsubsection 4.1.4. 
The energy density scalar behaves as W^(t) — VFo '^e(°ja2|^°'* • 

Second special example: if(a)(to) = 

When Ha{to) = 0, one obtains the simplified formulas 



Hia)it) 
<t) 

a(to) 



— tanh 
3 



^(l+/i||)(t-to) 



cosh 



-(H-/X||)(t-to) 



KX{t) 



A 



2yU|| — l) cosh 



.-2 



-(l+/i||)(t-to) 



The equation for Y can be transformed into the Legendre equation 

A2 



1 - z')Y"{z) - 1zY\z) + Y 



1-^2 



0, 



where 

z = tanh 



[l+fl\\){t-to) 



viv^l) 



5(2/i||-l) 
4(l+/x||)^ 



4(l+/i|| 



Thus, one obtains 
c(t) 



cito 



\-z\t) 



4{1+M||) 



[C{Pl{z{f)) ^ C^Ql{: 



(73) 

(74) 
(75) 

(76) 

(77) 

(78) 



where V^{z) and Qt{z) are the associated Legendre functions of the first and second 
kinds, respectively (see, |36], 8.1.1). Ci and C2 are the constants of integration 



Wr{0) 



(QtnO) - JQtiO) 



Co 



Wr{0) 



JVM - (^')'(O) 



(79) 
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where 



wr(Q) ^ vmm\^)-ivt)'{mm 



2^^r(|z/ + |A + i)r(|z/ + |A + |: 



is the Wronsky determinant at 2; = (t = to)) r(g) is Gamma-function, 



J 



(80) 
(81) 



and 



P^^(O) =2V-5^'2 

(n')'(o) 



r(iz.+ iA + i) 



r(|z/-iA + i) 



COS 



7r(z/ + A) 



r(iz/+iA + l) . 

TT 2-—^ f- ^ sm 



r(|^-|A + i; 



vr(z/ + A) 



Q'(0) 



.A-i_ir(i/^ + iA + i) 



7r2 



sm 



TT 



(Z. + A) 



(Q^)'(0)^2V^(^^ + ^^ + ^) 



cos 



7r(z/ + A) 



(82) 



Third special example: A = 
The solution of the problem is 



KXit) = (1 - 2^11) Hl^ito) X-^^l+'^N) , H^a){x) = ±H^a){to) 



X 



(83) 



^ = ^^"'^"W, H^a){t) = ^^^, ^(t)^l + /7(„)(to)(l + /X|i)(t-to). (84) 
The equation fl55l) reduces to the Euler equation 



z'Y"{z) + Y 



and c{t) reads 



where 



c(t) = c(to)2: 



0" 



/i^ + 12/i|| -4 
2(1+ /ill) 



2{1 + M||) 



+ 



5(l-2/i||) 

4(l+/i||)2 

c(to) 



/I 
2 



-4 



(85) 



(86) 



7) 



This solution may not appear in the model without magnetization, since the nec- 
essary condition 2H'^^^(to){2fi\^ ^ 1) + i^Bq^Mq + Bq) = is not valid when Mq = 
and = 1. Note that the asymptotics for a{t) and c(t). 



/'|| + ^M^|+12m||-4 



a(t ^ 00) oc t'+'^ii , c(t^oo)oct '''+'^11^ , (88) 

coincide only in the critical regime — | or in other words, when ^ \. Thus, 
when A = the Universe does not isotropize for arbitrary 7^ i. 
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5.2 Second example of cosmological evolution: longitudinal 
quasi- vacuum 

When W + P\\ = 0, the right-hand sides of (l28l) and of (l30l) coincide, and their 
left-hand-sides give 

ca a , ^, 

— = -. 89 

ca a 

One should distinguish two cases: first, when d(t) ^ and H(^a)(to) 7^ 0, second, 
when d = 0. 



5.2.1 a(t) = const = a(to) 

The Einstein equations, supplemented by the equations of state P\\ = —W and 
Ptr = (7 — formally admit the solution a{t) = a(to), if, first, the quantities 

P||, Ptr, W and X are constant and are hnked by the relation A + k{W + X)) = 0, 
second, the equation for c{t) has the form 

c + c{t)Q = 0, Q = k{2- X)X{a{to)) --fA = const. (90) 

When Q > 0, c{t) oscillates harmonically with the frequency y/Q, and this model is 
singular. When Q < 0, c{t) behaves exponentially. At Q = c{t) is a linear function 
of time and this model can be effectively reduced to the Minkowski spacetime. 



5.2.2 ay^O, Ha{to) ^ 
In this case the consequence 



gives c(t) readily in terms of a{t) and Hia){t) 



c{t) = c{to) 



Hia){t) ait) 
H(a)ito) a{to) 



(91) 



Asymptotic behaviour of c{t) at t — > 00 is the same as for a{t) (i.e., the Universe 
isotropizes) , when H(^a){t) — >■ const. The optimal strategy to obtain the solution is 
now the following. First, we find W from the conservation law fl3Tl) transformed into 



A (,^W(x)) = -x^^-^^, (x^^(x)) . 



d 



dx 

Second, we solve the equation for H'^{x 

d 



dx 



X 



— i^x'Hl))=A + K[Wix) + Xix)] 



(92) 



(93) 



which is the direct consequence of the equation (l3Up . Third, we consider the solution 
of fl5^ for a{t), and then return to the solution fl9T]) for c{t). 
When 2(7 - 1)^ ^ 1, the solution of ([92]) is 



W{x) = WqX-^^ + 



19 



1] 



-27 
X ' — X 



(94) 



When 2(7 — 1)^ = 1, W{x) behaves according to the formula 



W{x) 



X ^ 



Wn 



2^ 



■ log X 



In the first case, when 2(7 — 1)^ 7^ 1, the solution of (|93l) takes the form 

A 



+ x-^Li + x-^L2 + x-^^Lg + x'^-iL. 



'4 5 



(95) 



(96) 



where 



Li — Hf^^^ito) — — — L2 — L3 



1 



kBIKi 



kB? 



L 



4 , 



2(3 - 27) 



2Wq 



+ K2 



[2^(7 - 1) - 1] 



nB^nK2- 



^°^^'(e-l)[2^(7-l)-l]- ^^^^ 

Analogously to the case of paramagnetic / diamagnetic dust, we attract attention 
to the fact that generally H(^a){x) is not a monotonic function and has several ex- 
tremums. This means that in the expansion of the Universe there are time intervals 
characterized by acceleration and deceleration. Asymptotic behaviour of this func- 
tion is predetermined by the parameters A, 7 and ^. The novelty in comparison with 
the previous analysis is that a new parameter 7 is involved. To obtain a{t) we have 
to solve the equation fl5^ with H(^a){x) given by fl96l) . To illustrate our conclusions 
let us consider two particular cases. 



5.2.3 Slow relaxation 

When ^ ^ 00 one obtains a model in which relaxation time for the magnetization 
is much larger than the characteristic time of the Universe evolution. Consider for 
simplicity Li = and L3 = 0. It is possible if 



Mn 



2Wo(l-7) 



A 



1 



Ha){to) = i, + ^^Wo{2-^)--KB',. 



Bo ' -w— 3 
With such a choice of the initial parameters we have 



[a) 



X 



l!^ + ^Wo{2-j)x-^-^KBix-\ 



(a) 



kB?:x-^ I 1 



X 



(98) 

(99) 

(100) 



Here the definition xl = ,^^^2-7) used, and the assumptions 7 < 2 and B'^ > 
1^0(2—7) are made. At the point x=x^ the function if(a)(x) reaches its maximum 
with 



A KBi 
^3 + 2^' 



(101) 
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and tends to the value a/4 asymptotically at t — oo. Integration of yields 



' a{t) 



cosh 



2\/^ (t-to) 



'-^^(a)(io)sinh 



+ 



A 



(102) 



Thus, such a model gives a{t), -ff(a)(^) and c(t) in terms of elementary (hyperbolic) 
functions and is convenient for qualitative analysis. 



5.2.4 Transversal stiff matter, 7 = 2 

When 7 = 2, P(tr) = W, i.e., the matter behaves as a stiff one. It follows from (|92|) 
that 

Wix) + A'(x) = a;-^(iyo + ^o) , = ^So(Mo + So) , (103) 
and (!30|) yields 

^(a)(a;) = 



^(a)(^o)-| + «:(W^o + Xo) 



x-4fi:(l^o + Xo) . (104) 



This means that the parameter ^ becomes hidden if we consider the functions iJ(a)(t), 
a{t) and c{t). Nevertheless, it appears when we calculate W{x): 



W{x) = Wox-^ + ]^BlK2 



^-^ _ a;-(2+i 



(105) 



To illustrate the problem arising in a particular case if(a)(to) = mentioned in the 
beginning of Subsection 5.2, consider now the solution fll04p with a special choice 
of initial data: iJ(a)(to) = and h{Wq + A'q) = ^. Then one obtains 



a{t) = a(to) cosha 



(t-to) 



tanh 



(t-to) 



(106) 



and 



c(t) = const cosh 







sinh 


2/f (t-to) 



(107) 



The model is self-consistent, if c(to) = 0. The const can be chosen from the 
isotropization condition a{t 00) = c{t 00), i.e., const = a{to). 
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6 Constant relaxation parameter 
6.1 Evolution of magnetization 

When r takes a constant value, tq, we can readily solve ( l27ll in quadratures to get 



M(t) = exp 



t - tc 
To 



M{to) + 



— -1 



* dt ft -to^ 

exp 



to a'^{t') 



To 



(108) 

To obtain some analytical results consider the evolution equations with a{t) = 
a{tQ)e^''^^~^°\ where Hq is a constant. Such a two-dimensional de Sitter-type ex- 
pansion is possible, e.g., when 



P\^it) = X{t), A = 3H^ 



(109) 



In this case the total longitudinal pressure, V\\ = P\\{t) — X{t), vanishes. When 
2HoTq 1 the formula fllOSp for the magnetization and formula fl32|) for the X[t) 
read 



M{t) = exp 



t-tp 

To 



M(tp) + 



Bp 



2HpTp) 



exp [-2Hpit - to)] 



Bp 



(1 - 2HpTp) 



and 



+ 



Mjtp) 
Bp 



52 



1 + 



1 



{2HpTp 



{2HpTp - 1 



exp [-4iJo(t - to)] 



(110) 



exp 



(2Hp + tp) 



'111) 



respectively. The asymptotic behaviour of X(t) is dominated by the exponential 
function exp [—4:Hp(t — tp)] when > 2Hp, (i.e., the double relaxation time is 
less than the characteristic rate of expansion I /Hp ), and by another exponent 
— (2Hp + {t — tp) , when ^ < 2Hp. Note that a special case with X{t) = 



exp 

exists for 

M{tp) = -Bp , 2/7o/i||ro = 1 • 

The magnetization M{t) for such decreases exponentially 

M{t) = -So exp [-2Hp{t - tp)] . 
In the resonance case, when — = 2Hp, the formulas fllOSp and 



give 



M{t) = exp 



X{t) 



Bl 



exp 



Tp 



t-tp 
Tp 



■{t-tp) 



M{tp)-^{l--\ {t-tp 



Tp 



\ I M{tp) 
Bp 



t - tc 

Tp 



(112) 
(113) 

(114) 
(115) 
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Similarly to the case of time dependent relaxation parameter, the function X{t) can 
be monotonic or non-monotonic, depending on the value of the parameters M{to), 
Hq, tq. To illustrate this fact consider the resonance case = 2Hq, when /i|| > 1. 
If the following inequality takes place: 



5, 



^(^o)>^(l-3/i|i), 
2/^11 



(116) 



the function X{t) decreases, passes through its zero value, reaches a minimum at 
the point 



with 



4/i|| 



exp 



(3/i|| -l)Eo + 2/X||M(to) 
- l)Bo 

(3/i|l -l)5o + 2/i||M(to) 



- l)Bo 



< 



(117) 



1181 



and then increases and tends to zero asymptotically. For the diamagnetic medium 
we have to change the signs of inequality in flll6p and in (11181) . i.e., the function 
X{t) reaches the maximum. 



6.2 Third example of cosmological evolution: hidden induc- 
tion 

Consider the special case when X{t) = despite the magnetic field is no n- vanishing 
(see, also subsubsection 4.1.4). Assuming that P\\(t) = 0, we guarantee that the first 
equation in fllOQp is identically satisfied. Likewise, assume that P(tr) = (7 — 1)W. 
Then, Einstein's field equations (!28l) -( l30l) effectively reduce to the pair of equations 



2H. 



-W . 



K 



W + WHo{2-i + 1) + —rW = . 



The solution to the second equation is 

W{t) = H^(to)e-(2^+^)^°(*-*°) |l + 
thus. 



it) = c(to)e^°(*-*°) <^ 1 + 



2if2(27 + l) 



^(27+l)Ho(t-to) 



I _ g-(27+l)Ho(i-io) 



When t 00 one has the asymptotic relationship 



d 

dt 



log 



c(t) 
a{t) 



c a 
c a 



2H. 



g(27+l)//o(t-to) _ I 



(119) 



(120) 



:i21] 



[122) 



This means that the Universe isotropizes. 
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6.3 Fourth example of cosmological evolution 



Consider now the special case when V\\ = P\\{t) — X{t) = 0, a{t) = a(to) ^^p{Ho{t — 
to)}, A = 3Hq, and P(tr) = ^^P\\- For such a model the third Einstein equation 
( l30l) converts into identity, the second Einstein equation ( l29l) transforms into the 
equation for c(t) 



c + Hqc + c k{uj + l)X{t) - 2H^ 







(123) 



and the first one, (125]) . gives W(t) if c{t) is known. Note, that when = we have 
"transversal material dust", when = 1, the pressure of matter is isotropic. When 
= — 1, one obtains that V{tr) = -P(tr)+'^='P||=0. 

6.3.1 First special case uj = —1 

This case is the simplest, the solution of (1123^ is 



2c(to 



c(to 



Ho{t-to) , 

^3 



c(tc 



c(tc 



^2Ho{t-to) 



(124) 



The energy density can be found from the formula 
k{W + X) = -QHl \ 1 



2c(to) + ^ 



(125) 



where X is given by fillip . Initial value c(to) is connected with c(to), ^-nd Xq 
by the relation = Hq + which is a direct consequence of fl28|) . The 

Universe asymptotically isotropizes, i.e., a oc exp{Hot), c oc exp(ifo^)! and the total 
energy density W = W + X decreases as exp(— 3i/o^)• 
6.3.2 Second special case (cu ^ —1) 

For the special choice of the initial parameter M{to), which yields 



M(to) = 
the substitution 

with 



/i||(2/7oro-l)' ^' 2^||(2/7oro-l) 



rft:(^ + l)52(2/7oro/ii|-l)l^ 



A 



reduces the equation (11231) to 



(126) 

(127) 
(128) 

(129) 
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It is the Bessel equation (see, ^6], Eq. (9.1.1)). The solution can be expressed in 
terms of Bessel functions of the real argument when 



To > 



1 



2H. 



or To < 



2ifo/i|| ' 



for the paramagnetic medium, and when 



1 1 

To < TTTT or To > 



2Ho 



(130) 



(131) 



for the diamagnetic medium. This equation can also be reduced to the generalized 
Bessel equation for the imaginary argument iz (see, [36], Eq. (9.6.1)), when 



1 1 

< To < 



for > 1, and when 



1 



> To > 



2Ho 
1 

2^ 



(132) 



(133) 



2Hofi\\ 

for /i 1 1 < 1. In the last case c(t) can be expressed in terms of Bessel functions 
Ii,{z) = i^^Jyiiz). For simplicity we assume that u; + 1 > and tq is positive. Then 
the solution of (11291) is 



c(t) 



e 2 



Hoit-to) 



2Ho(t-to) 



(134) 



where Jiix) and J_j,{x) are the Bessel functions of the indices v = j and z/ = — |, 

4 4 4 4 

respectively. The constants Ci and C2 can be expressed in terms of initial data c(to) 
and cito): 



TT 



Co 



2^/2 

TT 



Ac{to)Aj\M) + (c(to) + ;^c(to) ) J_3(A) 



2^2 



4c(to)AJ4 {A) + ( c(to) + ;^c(to) ) J3(A) 



(135) 



(136) 



When t — i> 00 the argument of the Bessel functions in (11341) tends to zero, and we 
have the following asymptotic expression 



c{t 00) = C2 



2i 



,Ho{t-to) 



(137) 



where T is the Gamma-function. Thus, the de Sitter regime appears at t ^ 00. 
It is interesting that at t — ^ —00, when the argument of the Bessel functions tends 
to infinity, the corresponding formula 



c(t) 



V TT 



Cicos f Ae-2^"(*-*")- 



57r 



-C2 sin 



138) 
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demonstrates the fast quasi-harmonic oscillations of c(t) with the standard expo- 
nential damping. Note that M{t) vanishes when M(to) = and = 1. This case 
can also be described by the formulas ( 1134^ . ( 1135^ and ( 1136^ with 



A 



. 8^0^ J 



(139) 



and To 7^ jjj^. The energy density W as a solution of fl28l) inherits the dependence 
on time via the Bessel functions, we do not reproduce this expression here. 



6.3.3 Third special case 

When 2Hq^\\tq = 1, the term X{t) reads 

Mo 



X(t) 



1 + 



Bn 



(140) 



In this case the solution of the equation fll23p can also be represented in terms of 
Bessel functions 



c{t) = e 2 
where 



+M||){t-to) 



A 



k{uj + 1)5, 



1 + 



M(to) 
Bo 



:i4i) 



(142) 



,2i72(l + /i||)2 

Note that in the ferromagnetic phase fi\\ » 1 the argument of Bessel function in 
(11411) tends to zero much faster than in case (11341) . i.e., the isotropization in the 
ferromagnetic phase takes place faster. 



6.4 Fifth example of cosmological dynamics: 

non-homogeneous and non-hnear equation of state 

Consider now the special type of equation of state 

P(tr) = (7 - 1)W + . (143) 

Numerous non-homogeneous and non-linear equations of state of such kind are under 
discussion (see, e.g., [371 [38]). As before, we assume that A = 3Hl and a{t) = 
a{to)e^°^^~^'^\ For the equation of state fll43p the equation (J3T1) yields 

j^{W + + Hoi2j + 1){W + X) = -2HoXX , (144) 
whose solution reads 

W{t) = -Xit) + Li e-^»(27+l)(t-to) + g-4Ho(i-io) + ^3 ^-{2Ho + ^)(t-to) ^^45^ 
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Here the constants Li, L2 and L3 are given by 

Bl , 5oM(to) [i/oro(27 + 2A-l)-l] 



Li = W{to) + 



+ 



2 [iforo(27 - 1) 

i/oro/i||(27-l)-l 



1] 



/i||(27-3) [i7oro(27-l)-l] 
2ifoi"oAt|| - 1 



\Bl 



/i||(27-3) (2/foro-l^ 



HqXto 



+ 



L3- BoMito) [^^^^(27-1)-!] ' /X||(2iforo-l) [HoTo{2j - 1) - 1] 
Likewise, for c(t) we obtain 

c{t) ■ 



+ 



log 



KL2 
8HI 



Ho{t-to) 



1 _ g-4/^o(t-to) 



2^f(27 + T 
KL3 



1-e 



-/fo(27+l){t-to) 



1 _ g-(2J^O + -^)(i-io) 



(146) 
(147) 
(148) 



In the asymptotic regime at t ^ 00 one obtains c(t) c(oo)e'^°^*~*°'', where 



c(oo) = c(to) exp 



2i/n 



+ TT7- + 



i7o(27+l) 4/7o (2i/oro+i: 



(149) 



(150) 



Thus, the isotropization takes place, as it should. For t 
sup erexp onent ially. 



-00, c(t) decreases 



7 Discussion 

We have considered the simplest model of the one-dimensional relaxation of matter 
magnetization in a strong magnetic field in the framework of the extended Einstein- 
Maxwell theory applied to Bianchi-I cosmological model. We have shown that this 
model admits a set of exact analytical solutions, depending on the set of guiding 
parameters. Let us emphasize the main aspects of the obtained results. 

1. Analogy with extended (causal) thermodynamics. 

The key element of the extended Einstein - Maxwell theory in the context of 
anisotropic Bianchi-I model is the one- dimensional relaxation equation fl27|) for the 
magnetization M{t). The key element of the extended irreversible thermodynamics 
in the context of isotropic Friedmann model is a relaxation equation for the bulk 
viscous pressure a{t) (see, e.g., [H]). These two equations 

tM + M ={—-i\bJ and rd + a = -3C- (151) 

) \ a{t) J a 
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look similar. {( is a bulk viscosity coefficient). In both cases the rate of evolution 
of the scale factor a{t) predetermines the relaxation properties of M{t) or In 
both cases the function M{t) or a{t), appears in the right-hand-side of the Einstein 
equations as an element of the source term. In both models new degrees of freedom, 
activated in matter by the cosmological evolution, change the rate of expansion. In 
order to check this claim one can simply compare the results obtained for M{t) = 
and M(t) 7^ 0. As an example, let us compare the expressions (jH]), (p6l) . (1971) and 
(pTj) with those at Mq = 0, = 1, ^ = (i.e., at K2 = 0). The difference is 
that the magnetization adds a principally new second term in (JSj), describing the 
evolution of the energy density scalar, as well as the new last term in (l96ll . describing 
the rate of expansion in the cross-section x^Ox"^. Taking into account the formula 
fl9T|) . one can see that the modifications in H(^a){t) lead to the changes in the rate 

of evolution in the direction x^. Moreover, the presence of the new terms 
allows us to choose the phenomenological parameter ^ so that this term becomes 
of the leading order at x — >■ cxd in comparison with the terms x~^, and x~'^'^ in 
the formula (|96l) . In such a case just the magnetization predetermines the rate of 
cosmological evolution at t ^ 00, and the relaxation time r(t) = C,Ha{t) introduces 
a new expansion time scale. 

In both theories the relaxation time r is considered to be a function of cosmologi- 
cal time t and is a subject of modeling. In the extended irreversible thermodynamics 
the relaxation time is considered as r = -j^F, where ( = aW^ (in our definition of the 
energy density scalar) (see, e.g., [17] - [21]). When the function W{t) is obtained 
from the cosmological dynamics, the function T{W{t)) becomes an alternative rep- 
resentation of the function T{H{t)). 

The main difference of the results is connected with the fact that the Bianchi-I 
model is anisotropic, and M{t) is in fact a projection of the magnetization on the 
direction pointed by the magnetic field. In the isotropic Friedmann model a is a 
scalar describing the isotropic bulk viscous pressure. As a consequence, H{t) in 
the Friedmann model is positive and the right-hand-side of the relaxation equation 
for a is always negative. The sign of the right-hand-side of the relaxation equa- 
tion for M{t) depends on the sign of Bq, as well as on the sign of the difference 
(/i|| — 1). Respectively, the obtained magnetization may be positive or negative 
depending on the (random) initial value Mq, relaxation time and magnetic perme- 
ability. This option allows to consider a principally new situation, when magnetic 
field and magnetization are non-vanishing, nevertheless, the total magnetic source 
term X = \B{B + M) is equal to zero and disappears from the Einstein equations. 
Such solutions are discussed in the subsections 4.1.4. and 5.1.2. 

2. Monotonic and non-monotonic expansion 

Classical models with pure magnetic field are characterized by the non-negative 
source term X = ^B^ > 0, which decreases monotonically as a~^. The magnetiza- 
tion changes the situation: X may be positive, negative or equal to zero. Generally, 
X{t) is not monotonic function any longer, it may possess one, two or more ex- 
tremums. As a consequence of this behaviour, the function if(a)(t) is not monotonic, 
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thus, in the evolution of the Universe in the cross-section x^Ox"^ there are periods 
of (transversal) acceleration and deceleration. The simplest behaviour H(^a){t) is 
characterized by the presence of one minimum or maximum, and by the asymptotic 
de Sitter regime Ha{t) — >• (see, e.g., subsection 5.2.3). The behaviour of c(t) is 
also non-monotonic in this case. More complicated situation is characterizes by the 
solution for c(t), presented in terms of Bessel functions (see, subsubsection 6.3.2). 
The function c(t) behaves quasi-periodically, and one can expect that the number 
of periods of the longitudinal acceleration and deceleration is infinite. 

3. Guiding and resonance parameters 

The considered extended Einstein-Maxwell model is characterized by eight guiding 
parameters: ^ or tq, fi\\ — 1, Bq, Mq, Wq, 7, H(^a){to) and A. There are several 
underlined values of the parameter ^ (in the model of variable relaxation time) 
and of the parameter tq (in the model of constant relaxation time). The values 
^ = I and = 2Hq, respectively, are in fact resonance parameters, which appear 
in the integration of the differential equation for M{t). In such a resonance case the 
relaxation time r = ^H^J^{t) or Tq = jHq^ coincides with the characteristic time 
of the evolution of the function B{t) = F12 a^^{t), which provides the dynamics of 
magnetization. In case of resonance the function , as a part of M{t) (l34ll . has to 
be replaced by logx (see, (H^ ). It is very interesting to emphasize that in [21] the 
special value of the parameter q, q = ^, leads to the law r ~ H^^ for the relaxation 
of the bulk viscosity pressure. The value ^ = 1 is evidently the resonance value of 
the parameter ^, but it has another origin. When ^ = 1 the rate of change of the 

function H(a)(t) coincides with that of X(t). As a consequence, the function x~^^^'^^ 
in ( 1571) has to be replaced by x~'^loga; (see, ( l59l) ). The special values | = 2(7 — 1), 
^ = Ho{2-f - 1) (see, ([9lD)-([97D) and (fTl6D V ffT48D ) can also be considered as some 
analogs of the relation = 3ifo^^^i^ appeared in [i21j in the context of evolution 
of the bulk viscous pressure. The special value 7 = |, appearing in fl97|) . relates to 
the vanishing trace of the matter pressure tensor P1 + P2 + -P3 = 2-P(tr) + -P|| = 0. The 
special value = 2HqH\\ appears in the context of the vanishing X{t) (see, (11121) ). 
Finally, the model characterized by the constant values of P\\, -P(tr), W and X, 
discussed in the subsubsection 5.2.1 in context of the special condition Hi^a){io) = 0, 
also has an appropriate analog in |21j . 

4- Isotropization 

The models, in which the cosmological constant A is non-vanishing, isotropize at 
t —>■ 00. The exceptional case (see, (j88l) ) corresponds to A = 0. The first novelty 
of the obtained results is that at ^ = — | the magnetized matter can effectively 
redefine the cosmological constant (see, (1^ ). This case is exotic, since negative 
^ corresponds to magnetic instability and the magnetization increases with time. 
Nevertheless, when A = 0, this effect can in principle produce a non-vanishing 
effective cosmological constant. The second novelty is connected with the non- 
monotonic character of the isotropization. 
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8 Appendix: Variation of the tetrad vectors 

Let A^*^-) be the set of four tetrad four-vectors, whose index (a) runs over four 
values: (0), (1), (2), (3). Let A^^q-^ coincide with U\ the four-vector of velocity of the 
medium as a whole, and let A^^g^ coincide with A*, the director four- vector of the 
medium with uni-axial symmetry. The tetrad four-vectors are assumed to satisfy 
the orthogonality - normalization rules 

<?.,A(^,)Af,) = r^(„)(,) , r^W^'') Af„)A(«,) = g^" , (152) 

where r]{a){b) denotes the Minkowski matrix, diagonal (1, — 1, — 1, — 1). Since the 
tetrad four-vectors are linked by the relation containing the metric, we have to 

define the formula for the variation -j^- Varying the first and second relations 
( 11521) with respect to the metric, we obtain, respectively, 

Xk{b)SX^^^ + Afc(a)5A('^) = -A(„) A(^j,)(55(ifc , (153) 

SgP'^ = ^i^m [A?,)5Af,) + Af,)5A(^J . (154) 

The variation of arbitrary origin ^A^*^-, (not necessarily caused by the metric varia- 
tion) can be represented as a linear combination of the tetrad four-vectors: 

SXl) = XI)YJ'^ ■ (155) 

The tetrad tensor Y(^a) '^^ '^'^^ generally sjTiimetric. Using the convolution of fll54p 
with tetrad vectors, we obtain 

y(a)(6) ^ y(6)(a) _ ^gP^X^^'^Xf^ , (156) 

where we use standard rules for the indices, e.g., A^-'^) = rf^f^'^^^ Qq^X'^y Conse- 
quently, the symmetric part of the quantity Y'^°''>^^\ indicated as Z^°''^^^\ can be 
readily found: 

1 
2 

and the law (11551) reads now 



^(a)(b) _ -SgPIXl^'^Xf'^ , (157) 



5X1^ = ^Sg^'^ [a,(,)5; + A,(,)5;] + A(^,)Z(i)^) . (158) 
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Here is a skew-symmetric part ofY^^^, i.e., 2Z(^a){c) = ^(a)(c) — ^(c){a)- Therefore, 
the variation of the metric produces the variation of the tetrad, described by fll58p 
with vanishing skew-symmetric part Zi^a){c)- Thus, one finally has 

- - X,(,)5^ + X,(,)5; , (159) 



and we can use this formula for the variation of the four-velocity vector W = X^^^^ 
and for the variation of the space-like vector X* = X^'g-j . 
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